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Abstract 

Glasses are out-of-equilibrium systems aging under the crystallization threat. During ordinary 
glass formation, the atomic diffusion slows down rendering its experimental investigation impracti- 
cally long, to the extent that a timescale divergence is taken for granted by many. We circumvent 
here these limitations, taking advantage of a wide family of glasses rapidly obtained by physical 
vapor deposition directly into the solid state, endowed with different ’’ages” rivaling those reached 
by standard cooling and waiting for millennia. Isothermally probing the mechanical response of 
each of these glasses, we infer a correspondence with viscosity along the equilibrium line, up to 
exapoise values. We find a dependence of the elastic modulus on the glass age, which, traced back 
to temperature steepness index of the viscosity, tears down one of the cornerstones of several glass 
transition theories: the dynamical divergence. Critically, our results suggest that the conventional 
wisdom picture of a glass ceasing to flow at finite temperature could be wrong. 
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SIGNIFICANCE 


Does the glass cease to flow at some flnite temperature?” Answering this question -of 
pivotal importance for glass formation theories- would require ridiculously long observation 
times. We circumvent this infeasibility relating the (directly inaccessible) ultra-viscous flow 
of a liquid to the elastic properties of the corresponding glass, which we measure as function 
of its age. The older is the glass the lower is the temperature at which viscosity can be 
determined. Taking advantage of physical vapor deposition, we rapidly obtain a wide spec¬ 
trum of ages rivaling those of millenary ambers, enabling viscosity determinations at values 
as large as those pertaining to the asthenosphere. Our result ultimately rules out the flnite 
temperature divergence of the molecular diffusion timescale in a glass. 


INTRODUCTION 


A wide range of liquids can be cooled below the freezing point avoiding crystallization. 


Under such circumstances, viscosity increases dramatical 


falls out-of-equilibrium [iy^. Whether j^lO| or not 
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up to the point where the system 


16| this increase signals a divergent 


behaviour is a question hitting at the core of any glass 


an increasing number of experimental investigations 


ransition theory, and has attracted 


17l-l27|. 


The thermal history of a liquid engraves the properties of the resulting, amorphous mate¬ 
rial (glass), which continuously hikes down a rugged energy landscape to enhance its stability 


in search of lower minima 


28 


29| . From an opposite perspective, the way the equilibrium 


properties of a unique liquid are encoded in the peculiarities of different corresponding 
glasses, each endowed with different histories and stabilities, has always puzzled scientists 
and technologists js, 30|. The Active temperature Tf, for instance, unambiguously signals 
the out-of-equilibrium arrested state accompanying the glass transition {mI, and increases 
with the cooling speed, which is upper limited by the flnite thermal conductivity of the 
sample. Upon slowing down the quenching rate, the system gets trapped in lower potential 
energy landscape basins (inherent structures), but crystallization is also facilitated. In an 
attempt to rationalize such phenomenology, one way to define the glass transition temper¬ 
ature is as the Active temperature at a conventional cooling rate: Tg = Tf{—10 K/min). 
Correspondingly, the system falls out-of-equilibrium over a timescale of ~ 100 s, attaining 
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viscosities of 10^^ 10^^ poise. The impact of such an arbitrariness is mitigated by the 

observation that, for bulk materials, ordinary cooling rates are often confined in the range 
-0.1 to -10 K/min and Tf varies very little. Remarkably, this severely limits the free energy 
region accessible by the glass, at least on timescales useful for many potential applications. 


It has been recently shown 


28 


32 


34l |. however, that the limitations in exploring the bot¬ 


tom of the energy landscape can be circumvented synthesizing glasses directly well below 
Tg, avoiding the supercooling route by physical vapor deposition (PVD). These glasses are 
referred to as ultrastable, for their enhanced resistance against crystallization, and populate 


free energy regions as low as those of semi-fossilized resins 


17| aged for thousands of years. 


Indomethacin (IMC), selected for this study, is a prototype organic glass-former endowed 
with antipyretic properties, which, in its crystalline form, has poor water-solubility and 


bioavailability 


35| . limitations partially removed in the amorphous phase. Obtaining ul¬ 


trastable amorphous IMC on laboratory timescales does not improve water solubility, but 
enables a unique opportunity to approach the so called ’’ideal glass”, i.e. an amorphous 
phase with the same entropy of the undelying crystalline form. 

RESULTS 


Glass stability can be quantihed by determining Tf from the endothermic jump exhibited 


by the specific heat upon fast heating 


, due to the transformation of the glass into 


supercooled liquid (Materials and Methods). We present in Fig. 1 an assessment of the 
stabilities of IMC glasses prepared by PVD at six different substrate temperatures Tg^b in 
the range 190 K 290 K, as determined using differential scanning calorimetry (DSC). 

The Tj’s of the depositions testihes to their stability enhancement (compared to con¬ 
ventional glasses), which increases with the substrate temperature up to a maximum gap 
of 34 K with respect to Tg, achieved for the most ultrastable, lowest Tf case (Ts„b=266 K, 


in agreement with 


33|b The reason for such stability enhancement, attained by tuning 


Tguf,, has been traced back to the increase of the surface mobility, whic 


1 is believed to allow 


molecules to arrange themselves into lower energy configurations 


28 


32 


331. Specifically, for 


a given deposition rate, coming from the high temperature side, the (relatively large) surface 
mobility grants that the glass is ’’assembled” in equilibrium with its liquid, i.e. each layer of 
molecules is equilibrated before new molecules arrive in the deposition process. Within this 
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regime, the lower is the temperature, the lower the entropy of the system and the higher the 
stability. Below a given temperature the mobility is so slow that the structure is no longer 
in equilibrium with the liquid, and the system becomes increasingly unstable 

Since the shape of the inherent structures of the energy landscape defines the vibrational 
properties, a question arises on the possible connection between mechanical response and 
glass stability, in connection with recent experimental evidences in amorphous depositions 


and naturally aged glasses 


HE, 


38| . Specifically, previous studies |39j of IMC obtained by 


PVD reveal a sound velocity dependence on local substrate temperature which, in turn, has 
been shown to relate to the onset temperature of enthalpy release measured in calorimetry 
upscans Hi- To address this issue we first studied the acoustic properties of the most ul¬ 
trastable deposition by a combination of time domain pump-probe optical spectroscopy and 
frequency domain Inelastic X-ray Scattering (IXS). The former technique is a broadband 


version 


41 


42| of picosecond photo-acoustics (BPA) |^, which we recently developed to 


efficiently generate and detect longitudinal acoustic wave-packets in sub-pm thick transpar¬ 
ent layers. Specifically, by using white light pulses (450 650 nm), the energy dispersion 

and mode attenuation of the vibrational excitations have been simultaneously investigated 
in the frequency range 10 ^ 20 GHz. IXS, on the other hand, allows the study of the high 
frequency limit (THz), where the characteristic excitation wavelengths approach the average 
interparticle distance and the structural disorder becomes of crucial importance. Further 
details on both the techniques are given in Materials and Methods. The dispersion curve 
as determined by BPA in the low frequency regime, and by IXS in the mesoscopic regime, 
is reported in Fig. 2, along with the photoinduced differential reflectivity oscillations map 
and IXS spectra at fixed exchanged momenta (insets a and b), from which the dispersion is 
obtained. 

In the low-frequency region, the linear energy dispersion indicates the existence of prop¬ 
agating vibrational excitations corresponding to a single longitudinal acoustic-like phonon 
branch. Remarkably, the propagating nature of these excitations lingers up to the THz 
frequency regime. The presence of structural disorder, indeed, strongly attenuates acous¬ 
tic modes at the nanometer scale. Nevertheless, the existence of a well defined first sharp 
diffraction peak mimics a first pseudo-Brillouin zone and, though the vibrational eigenmodes 
are no longer plane waves, a well defined dominant wavevector still exists up to a zero-group 
velocity point corresponding to half of the pseudo-Brillouin zone boundary The ob- 


4 












tained sound velocity (the q = 0 derivative of the best sinusoidal £t to IXS data, dashed line 
in Fig. 2) is in excellent agreement with BPA data in the ten GHz range and with previous 
single (lower) frequency determination 371. 

The existence of a well defined dispersion at any lengthscales within a pseudo-Brillouin 
zone, suggests the potential sensitivity of the acoustic properties to the local topology of the 
inherent structures visited by the glass. Given its capability to access a sizable portion of the 
dispersion curve, determining sound velocity with superior accuracy (down to 0.1 %, about 
one order of magnitude better than ordinary Brillouin Light Scattering) BPA measurements 
have been extended to the other PVD samples in Fig. 3 (a and b). During the experiments 
the temperature was kept constant at T=295 K, i.e. well below Tg and above the T/ of the 
different depositions, to ensure that the glass state is arrested on the laboratory timescale. 

A clear inverse correlation between sound velocity and Tf arises. BPA experiments also 
enable to determine the hypersonic attenuation P with an accuracy of tenths of GHz, and, 
accordingly, determine the frequency dependence of the life-time of the vibrational excita¬ 
tions. As an example, the frequency dependence of P for Tf = 281 K and 301 K are reported 
(pink and blue, respectively) in Fig. 3 (c). Aiming to unravel any link with Tj, we integrated 
the attenuation over the explored frequency range and calculated the excess relative (and 
normalized) to the most ultrastable glass (Tf = 281 K), as reported in Fig. 3 (d). The 
relative attenuation shows, similarly to the sound velocity, a clear but opposite (direct) cor¬ 
relation with^. Acoustic attenuation in disordered materials is ruled by different physical 
mechanisms 42|. In the GHz regime, where a crystal-like picture of vibrational excitations 
holds, the attenuation is expected to be due to the anharmonicity of the interparticle interac¬ 
tion. The decrease of sound attenuation towards the lowest T/s discussed above, therefore, 
corroborates the idea that lower energy basins are more harmonic in nature. 


DISCUSSION 

The distinct thermo-mechanical correlation reported in this study has important impli¬ 
cations for the validity of those paradigms which can be uniquely benchmarked at the liquid 
side of the glass transition, since equilibration times dramatically increase below Tg prevent¬ 
ing direct explorations. Fragility, in particular, quantifies the steepness of the liquid viscosity 
(or equivalently the relaxation time, through Maxwell’s relation rj = Goot) at Tg, and it has 
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a pivotal role in controlling physical properties of the supercooled phase [2|, yj. Remarkably, 
the provocative idea that fast degrees of freedom characterizing glass dynamics well below 
Tg can be predictors for slow dynamical properties of the liquid state, such as the structural 


relaxation time and hence fragility, is rapidly gaining consensus 


45 


481. While density 


fluctuations completely decorrelate in a liquid, as a consequence of the ergodic sampling of 
different basins in the energy landscape, in a glass, a residual correlation exists, uniquely 
determined by the entire spectrum of vibrational eigenstates of the energy minimum where 
the system is trapped. This is the so called non ergodicity factor (NEF), the long time limit 
of the density-density autocorrelation function. Of interest here, it has been shown that the 
iquid fragility, m, can be determined by the low temperature behavior of the NEF, fq{T) 

, The NEF can be quantihed in different ways. Among them, by the relative sound 


46 


velocity jump occurring at the glass transition 49 


50| . namely: 


m 


= 7-(/-'(T,)-1)^140- 


( 1 ) 


where Cq and c^o are the liquid and glass sound velocity values, respectively, and 7 ~ 140 
expresses the above discussed correlation with fragility m verihed in ordinary IMG 


glass 


5l| . Critically, Cq is an equilibrium property of the liquid, while Coo is shown here to 


depend upon the very stability of the corresponding glass. This advocates the extension 
of the fragility concept to a temperature dependent steepness index jli| . I{T) such that 
I{Tg) = m, which we use here to obtain the equilibrium viscosity below Tg. The mechanical 
response of a glass of a given stability, indeed, determines the viscosity of its liquid at the 
corresponding Tf, via the dependence I{T = Tj) = I{c^{Tf)) in Eq. [1], reported in Fig. 4 
(a). 

Integrating this latter (see Eq. [TT] in Materials and Methods), an Arrhenius plot can 


ultimately be o 
IMG above T„ 


Dtained, shown in Fig. 4 (b). The super-Arrhenius behaviour documented in 


52| is reported as dashed line representing a Vogel-Fulcher-Tamman (VET) 


function. For T < Tg a. remarkable deviation is observed: the apparent fragility is larger than 
m, but falls below the VET expectation, and decreases with the glass stability, signifying 
a fragile-to-strong transition, quantihed in Fig. 4 (a). This result verihes the predictions 
of Kovacs and Adam-Gibbs-Vogel models from aged IMG 27(] and rationalizes very recent 
work in naturally long-time aged amber [w], which set an upper bound to the temperature 
dependence of the dynamics below Tg. It is also in line with the conclusion against the 
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VFT extrapolations in the glass [h, ^|, and relates to the observation of an additional non 
Arrhenius equilibration process in polymers 


54| . Moreover, the simultaneous decrease of 


both the acoustic attenuation and genera’ 


recent molecular dynamic simulations 


ized fragility reported here at low T^-’s, validates 
55| which put forward a direct correlation between 


kinetic fragility m and the degree of anharmonicity of the interparticle interaction potential. 

All together, these evidences syncretize on the scenario schematically depicted in Fig. 4, 
establishing a hrm link between the hypsometric characterization of the energy landscape and 
basin-specihc vibrational properties. Crucial to technological applications, the mechanical 
approach advanced here to assess stability is a non destructive one, as opposed to calometric 
determinations of Tf. Most important, we demonstrate that the state of the glass is totally 
identihed by Tj and T, and once the non ergodicity parameter is calculated by the sound 
velocity jump at Tj it is possible to determine the relaxation time from equation Eqs. ([1]) and 
oa. as shown in Fig. (|4]). The emerging protcol provides a unique way to capture essential 
features of a liquid attaining the structural arrest away from Tg, up to relaxation times 
beyond 10® s, corresponding to exapoise viscosities. In conclusion, we circumvented here 
the major hindrances to probe equilibrium properties below the ordinary glass transition 
temperature: crystallization tendency and unrealistic observation times. The isothermal 
measurement of mechanical properties for glassy phases with different Tf reveals the increase 
of the longitudinal elastic modulus and decrease of the acoustic loss below Tg, which we 
connect to a breakdown of the VFT law. This may lead to reconsider old paradigms based 
on dynamical divergency to ever reach a dehnitive understanding of the glass transition. 


MATERIALS AND METHODS 

Sample growth and thermal characterization and treatment 

Indomethacin (IMC, 99% purity, Tg= 315 K and T^ (7 form)= 428 K) crystalline powders 
were purchased from Sigma-Aldrich. Twin hlms of ultrastable IMC glasses were grown by 
vapor deposition on Si (100) substrates and on differential scanning calorimetry A1 pans at 
deposition temperatures spanning from 190-to-290 K. An effusion cell hlled with IMC was 
heated to achieve the desired deposition rate of 0.1 nm/s, as measured by a quartz crystal 
microbalance. When this rate was attained, the shutter was removed to start deposition. 


7 










The thickness of the hlms ranges from 1.5-2 /rm for the BPA experiments to 50 /rm for the 
IXS measurements. All samples were stored in vacuum-sealed bags with desiccant in the 
freezer to minimize aging prior to the acoustic measurements, which were carried out several 
days after the deposition. Differential scanning calorimetry was performed immediately after 
removing the samples from the growth chamber 5^. A DSC Perkin Elmer 7 was used to 
monitor the power absorbed/released during heating scans at a rate of 10 K/min on IMG 
hlms of 2-10 mg. The hrst scan typically corresponds to an ultrastable glass, while the 
second one is characteristic of a glass obtained by cooling the liquid at -10 K/min. The 
hctive temperature is measured, evaluating the intersection of the enthalpy curves of the 
PVD glasses, i.e. the integral of the heat capacity upscans, with the supercooled liquid line 


extrapolated according to a quadratic ht 


32| . Notably, the extrapolated value of the hctive 


temperature does not depend on the heating rate and directly rehects the structural state of 
the glass. The lower values of the hctive temperature of ultrastable glasses signify a smaller 
enthalpy content compared to the glasses cooled at ordinary rates. 


Inelastic x-ray scattering 

The IXS experiments were performed at the beam line ID28 of the European Synchrotron 
Radiation Facility (ESRF). The experimental observable is the frequency spectrum of the 
scattered X-ray intensity, proportional to the dynamic structure factor S{q,fkj), where q 
is the exchanged momentum dehned by the scattering angle 6 and by the wavevector of 
the incident photons ki as q = 2kiSin{9/2). A 8-analyzers bench was utilised, operating in 
horizontal scattering geometry. The sample was mounted with the substrate parallel to the 
scattering plane, tightly focusing the beam down to 20/im in the vertical direction to avoid 
parasitic scattering from the substrate. With such a geometry the probed phonon direction 
is orthogonal to the growth. Energy scans were performed at constant q values in the range 
1 to 8 ~ nm“^, corresponding to half of the pseudo Brillouin Zone (q=14.7 nm~^). The 
q resolution was determined by slits placed in front of the analyzer was set to 0.25 nm~^. 
The scanned energy range was —30 < hw < 30 meV, where hw = Eq — E is the energy 
transfer, with Eq and E being the energy of the incident (23.725 eV) and the scattered 
x-ray photon; each scan took approximately 480 min. Using the (12 12 12) reflection for 
the Si monochromator and crystal analyzers the overall energy resolution was 1.5 meV the 




FWHM of instrumental function (black line under the spectra in Fig. 2). Measurements 
were performed at room temperature (T= 295 K) on a sample grown by PVD at substrate 
temperature T=266 K and deposition rate of 0.1 nm/s, 50pm thick. The dynamic structure 
factor, S{q, Itw), contains information about the sound dispersion fl{q) and attenuation F(g) 
which can be extracted by a damped harmonic oscillator model (DHO): 


S{q,uj) = S{q) 




(1 - fq) 


n{qfT{q) 


TT 


hjJj/ikBT) 


( 2 ) 


2 — 

were the Bose factor accounts for the quantum nature of the probed excitations. The two 
terms in Eq. [2] represent the elastic and the inelastic components of the spectra, respec¬ 
tively. Eq. [2] was htted to the experimental data after a correction for the hnite instrument 
resolution [57 1. 

The IXS sound velocity measurement of the ordinary glass reported in Fig. 3 (b) accounts 


for a 0.8% positive dispersion occurring in the THz regime 51| 


Picosecond photoacoustics 


Time-domain measurements of longitudinal sound velocity and acoustic damping in the 
hydrodynamics limit were performed by broadband picosecond photoacoustics. This pump 
and probe technique is based on the generation and detection of coherent vibrational exci¬ 
tations by means of ultrashort laser pulses. The setup is built on a regeneratively amplihed 
Ti:sapphire laser producing 50 fs, 4 mJ pulses at 800 nm with 1 kHz repetition rate. The 
output is split to generate both the pump and the probe beams; the former, after passing 
through a delay line, is focused onto the sample on a ~ 100 pm spot diameter with energies 
up to 5 pJ. The PVD glasses of IMG, 1.5 pm 2 pm thick, have been deposited onto 
a crystalline Si substrate at six different substrate temperatures (T= 190, 210, 236, 266, 
285, 290) and coated by a 15 nm thick layer of Ni. As the pump impinges on the sample, 
it is partially absorbed by the metallic coating, which thermally expands launching very 
short strain pulses, i.e. longitudinal acoustic wave packets in the underlying glass, with a 
characteristic spectrum extending from few to hundreds of GHz 431. 

The photo-generated acoustic wave packet travels inside the sample along the growth 
direction. A portion of the probe pulse is reflected at the metallic surface, while the trans- 
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mitted component interferes with the light scattered from the traveling density fluctuation 
associated with the acoustic wave packet. By monitoring transient differential reflectivity 
(with and without the pump pulse) as function of the time pump-probe time delay, 

information can be gained on the phonon spectrum. Specifically, at any given probe wave¬ 
length A, the phonon frequency is set by the Bragg condition corresponding to a stimulated 
Brillouin scattering event: 



where v is the longitudinal sound velocity, q the exchanged momentum, n{X) the index of 
refraction of the glass, and is the incident angle with respect to the metallic coating. In 
our realization, the probe pulse consisted of a broadband white light continuum, generated 
from 2 mm thick Sapphire crystal plate and hltered to select wavelengths in the range 450 
650 nm. The broadband probe pulse’s reflection was dispersed by a 150 g/mm grating 
and monitored by a CCD, such that a number of two color pump and probe experiments 
were simultaneously recorded. Transient reflectivity data were first reduced by subtraction 
of an exponential thermal background. Each channel achieves a sensitivity up to ^ Ri 10“®. 
In order to get rid of an independent determination of the refractive index n, required in 
Eq. [3l the measurements were repeated at different scattering angles (3. A selection of the 
reflectivity oscillations as function of time and probe wavelength is reported in color map 
in Fig. 2 together with the signals at selected wavelengths. The signal was modeled with 
the damped harmonic oscillator in time domain, whose frequency can be estimated with an 
accuracy of 0.1 %. The sound velocity was determined from the linear fit of the modes’ 
frequency v as function of the exchanged momentum q over all the broadband probed range. 
The acoustic attenuation, indeed, was directly extracted from the damping of the reflectivity 
oscillations with 6 % accuracy. 


Temperature dependence of the relaxation time 


The so called kinetic fragility quantihes the temperature behaviour of the relaxation time 
r(T) at a conventional hctive temperature {Tg) where r =100 s. Namely: 

dlogiQT 


m = 


diTg/T) 


( 4 ) 


T=Ta 
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In order to describe the low temperature dependence of the non ergodicity factor fq(T) in 
the low q limit, a dimensionless steepness index can be introduced: 

Of 


a = — lim 


1^0 d4r 


(5) 


T=0 


Remarkably, a exhibits a strong direct correlation with the kinetic fragility m, as it has 


been demonstrated experimentally |46 


£3, l58|. 


59| and by numerical simulations ^ for a 


sizable number of glass-formers. Accordingly, the faster is the T dependence of fq{T), i.e. 
the higher is the value of a, the more fragile is the glass-former. Below fragilities as high as 
one hundred, the correlation is linear: 


m = ja 

with 7 140. In the harmonic approximation, the T-dependence of fq{T) reads: 

1 


f,{T) = 


1 -|- ex- 


( 6 ) 


(7) 


The value of fq(Tg) can be conveniently determined from the glass (coo) —>■ liquid (cq) sound 
velocity jump around the glass transition, through the expression: 


MT.) = 1 - 


By evaluating Eq. Hat T = Tg, and combining it with Eq. IH it immediately follows: 


m = 7 — 




C2 - 

^OO 

*^0 


( 8 ) 


(9) 


Beside the measurement temperature, however, the glass velocity also depends upon stabil¬ 
ity, i.e. on the hetive temperature Tf. This can be taken into account by generalising Eq. [9] 
as: 


dlogioT 


d{Tf/T) 


= 7 


cl{T) 


T=Tt 


= I{Tf) 


( 10 ) 


T=Tf 


cliTf,T)-cliT) 

where the notation Coo{Tf,T) explicitly indicates the above mentioned dependencies on 
Active and measurement temperatures, respectively. It is worth mentioning that the T- 
dependence, in Eq. [TU] (i.e. on the very temperature where the sound velocity jump is 
evaluated) is negligible compared to the Tj dependence brought about by the glass sound 
velocity. Once integrated in 1/T, Eq. [TOl provides the relaxation time below Tg through: 


log,oT{l/Tf) = %ior(l/T,) + Tf / /(T')d(l/T0 


( 11 ) 
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with co(Tg) = 1354 m/s estimated by Brillouin Light Scattering [6^, and the c(T) depen¬ 


dence is from 


37l |. In Fig. 5 we sketch the evolution of the sound velocity as function of 


temperature across the glass transition, for two glasses of different hctive temperature. 
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FIG. 1: Fictive temperature of PVD glasses as function of the substrate temperature during 
deposition, determined from DSC upscans, directly relates to the enthalpy content, hence qnanti- 
fying the glass stability (see Materials and Methods for more details). Typical values obtained for 
conventional glasses for two different cooling rates are also indicated (dashed lines) 

indomethacin probed by brillouin light scattering. Phys. Chem. Chem. Phys. 16, 14206-14211 
(2014). URL http://dx.doi.org/10.1039/C4CP01427H. 
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FIG. 2: Longitudinal acoustic phonon dispersion in ultrastable PVD glass {Tguh=2GQ K), investi¬ 
gated by IXS (A, dashed line is the sinusoidal best fit) in the THz, b y B PA (o) in the GHz regime, 


along with a previous Brillouin measurement (□) at lower frequency 


371 ]. (a) Transient reflectivity 


BPA data, each Xprobe couples to a phonon of well defined wavevector, visualized as an oscillatory 
signal in time domain, (b) IXS spectra, acoustic phonons for selected wavevectors, Q, probed in 
frequency domain. 
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FIG. 3: Sound velocity and damping for different stabilities, (a) GHz-dispersion curves of the 
longitudinal acoustic phonons from isothermal BPA measurements (T = 295 K). (b) Sound ve¬ 
locity (linear fit to the dispersion curves), exhibits a correlation with Tj, showing a maximum for 
Tf = 281 K. Error bars smaller than the symbol. The values for conventional glasses (Tf = 307 K 


and Tf = 315 K) are from 


37 | and 


5l[ |. respectively, (c) Acoustic attenuation, F, for depositions 


with high (Tsub = 266 K, pink) and low (Tsub = 190 K, blue) stabilities as function of the phonon 
frequency, (d) Normalized excess acoustic attenuation relative to the most ultrastable glass, ob¬ 
tained integrating over the explored frequency range and normalizing to the Tf = 281 K case. An 
inverse correlation respect to that in sound velocity is observed. 
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FIG. 4: Viscous flow of a liquid and stability of its glasses, (a) Apparent fragility I{T), determined 
from the sound velocity jump in glasses with different Tf, indicates a fragile to strong transition 
below Tg. Conventional fragility m = I{Tg) is also reported, (b) Arrhenius plot for T < Tg, 
obtained integrating data from panel (a). A correlation is found between the depth of the inherent 
structures in the energy landscape (grey line) probed by glasses of different Tf and the mechanical 
response (sketched by time domain damped oscillations). The VFT behavionr obtained from above 
Tg data [5^ is also shown (dashed line). 
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FIG. 5: Temperature dependence of the sound velocity for two glasses of different Active temper¬ 
ature, Tf, along with the unique liquid and supercooled liquid values. The shaded line indicates 
the region of the highly viscous supercooled liquid, where direct observations are prevented by 
impractically long equilibration timescales. 
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